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Preliminary Structural Sizing of a Mach 3.0 
High-Speed Civil Transport Model 

1.0 Introduction 

The development of future high-speed civil transport aircraft 
conceptual and preliminary designs requires the integration of 
inputs from all of the traditional aeronautics disciplines. An effort 
currently underway at NASA Langley Research Center to integrate 
these disciplines is focusing on the conceptual design of a High- 
Speed Civil Transport (HSCT) to provide a testbed for the 
development of the multidisciplinary methodology needed for future 
vehicle design. The success of this venture is dependent on the data 
management system used by the design system and its ability to 
interface with the various discipline-oriented application programs 
for the exchange of data in a fast and accurate manner. 

Obviously, the development of a multidisciplinary, integrated 
design system is an immense task when considered overall, but the 
task is somewhat more manageable and the data requirements are 
more accurately defined when approached from a discipline level 
which acknowledges the inputs and outputs of the associated 
disciplines. These discipline oriented studies could be very 
informative when used in defining the application programs for the 
integrated system and the definition of the data base to support the 
input/output data requirements. In this regard, an analysis has been 
performed pertaining to the structural resizing of a candidate Mach 

3.0 HSCT conceptual design using a computer program called EZDESIT 
(Ref. 1). 

EZDESIT is a computer program which integrates the PATRAN 
finite element modeling program (Ref. 2) to the COMET finite 
element analysis program (Ref. 3) for the purpose of calculating 
element sizes or cross sectional dimensions. Bar, beam, and plate 
type structural elements are sized to support calculated stress 
resultants. The element sizing criteria consist of yield and ultimate 
strength, minimum gage, and local buckling constraints. The 
program is written in FORTRAN and operates on a CONVEX computer 
under a UNIX operating system with an interactive user interface. It 
is menu driven with interactive user queries which allows the user 
to translate PATRAN output to COMET input, structurally size the 
finite elements, and then translate the COMET output to a PATRAN 
postprocessing input. User instructions for menu selections and 
response to the interactive queries are given in Reference 1. 
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The purpose of the present report is to document the procedure 
used in accomplishing the preliminary structural sizing of a Mach 

3.0 HSCT model and present the corresponding results. In describing 
the procedure, the data requirements are discussed in detail in order 
that they may be evaluated with respect to the data management 
goals of the vehicle design project. Furthermore, those modeling 
requirements unique to EZDESIT are documented as they relate to the 
creation of the finite element model by PATRAN (Ref. 2) since the 
PATRAN neutral file is required input for EZDESIT. 

2.0 Foundation of Configuration and Aerodynamic Loads 

The external configuration of the baseline Mach 3.0 HSCT 
aerospace vehicle concept used in the present investigation is fully 
described in Reference 4. This concept, shown in Figure 1, is a 
blended wing-body configuration with a modified platypus nose, a 
highly swept inboard wing panel, a moderately swept outboard wing 
panel, and curved wingtips. This configuration was primarily 
selected for structural modeling studies because it satisfies the 
aerodynamic and payload requirements necessary for efficient 
operation within the civilian air transportation industry. During 
this early phase of the conceptual design only minor consideration is 
usually given to structural problems (i.e., depth limitations due to 
slenderness and/or contouring of the external shape). 

The aerodynamic loads are believed to be consistent with the 
load cases used for the study described in Reference 4 along with 
the associated baseline weights and performance data. These data 
are included in Table A-1 of Appendix A. A program called WINGDES 
(WING DESign) (Ref. 5) was used to compute the aerodynamic 
pressures for each of the load cases. Due to the limitations of 
WINGDES, only pressure distributions on the wing were considered in 
the analysis. 

3.0 Description of Finite Element Model 

A structural arrangement was generated within the confines of 
the external shape defined by Table IV of Reference 4 and, for the 
sake of completeness, is presented in Table B-1 of Appendix B. A 
modified version of the structural arrangement for the wing and 
fuselage configuration is shown in Figure 2. A brief discussion of 
the routine for creating this structural arrangement and its 
application to the generation of a wing structure is given in 
Reference 6. This arrangement was used to develop the PATRAN 
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Phase 1 geometry (grid points, lines, and patches) that is required 
for the generation of PATRAN Phase 2 geometry describing the nodes 
and elements of the finite element model. For this model, the 
network of lines and patches were defined by a set of grid points 
that occur at the intersection of the ribs and spars in the wing and 
the frames and longerons in the fuselage. Furthermore, the nodes of 
the finite element model were chosen to be coincident with the grid 
points of the geometrical model and, therefore, the representation 
of the finite elements is identical to the structural arrangement as 
shown in Figure 2. 

(a) Structural Model 

The structural elements of the finite element model consist of 
577 rod, 229 beam, and 508 plate elements that are interconnected 
at their end or corner points by 398 nodes. The rod elements 
represent the caps and web stiffeners for the ribs and spars and the 
longerons of the fuselage. The beam elements represent the 
fuselage frames. The 508 plate elements are a combination of 32 
triangular and 476 quadrilateral plates which are capable of 
representing both the membrane and bending stiffnesses of the cover 
panels for the wing and fuselage and the shear webs of the wing's 
ribs and spars. The triangular elements are required in modeling the 
rib structure at the leading and trailing edges of the wing and in 
some transition areas (i.e., wing-fuselage intersection and wing 
crank). The wing carry-through structure is a projection of the root 
rib and intersecting spars of the wing box through the fuselage in a 
direction normal to the airplane's plane of symmetry and is idealized 
accordingly. For convenience, the initial property identification 
numbers for each set of elements were assigned during this phase of 
model generation. The finite element model had a total of 2163 
degrees of freedom. 

Unique EZDESIT input requirements are satisfied during the 
formulation of the elements of the model; namely, assignment of a 
name or names to one or more structural components (i.e., wing 
upper surface, ribs, fuselage, etc.) and element numbering in the 
required sequence. EZDESIT requires that at least one component be 
assigned a name (Ref. 2, page 15-8) but good modeling procedures 
suggest naming all components in detail. For example, individual 
names are given to the wing's upper and lower surfaces, shear webs 
of the ribs and spars, and caps for the ribs and spars as opposed to 
giving a single name to the entire wing assemblage. This naming 
practice is particularly advantageous when assigning loads, 
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constraints, and material and section properties. Furthermore, the 
elements contained in the PATRAN neutral file must be numbered 
consecutively beginning with the rods and continuing with the 
beams, triangular plates, and quadrilateral plates, respectively. If 
necessary, the elements can be renumbered using the PATRAN 
"RENUMBER" command given on page 5-28 of Reference 2. The 
number of each of these elements, which is required user input for 
the EZDESIT execution, is obtained using the PATRAN "SHOW, MAX2" 
command. 

(b) Model Constraints 

For computational efficiency, only half of the HSCT was 
modeled and boundary constraints were applied to the nodes in the 
plane of symmetry to represent the remaining half. The rigid body 
motions were restrained by applying the necessary constraints to 
those nodes of the wing carry-through structure that lie in the 
airplane's plane of symmetry (namely; constraints in the Z-direction 
at two nodes along the X-axis to prevent vertical and pitching 
motion and a single constraint in the X-direction to prevent axial 

motion). 

(c) Material and Section Properties 

This initial study is based on two finite element models 
composed entirely of 2024-T81 aluminum or Ti-6AI-4V titanium, 
respectively. The resizing routine used in this study is presently 
restricted to structural components made of isotropic materials. 
However, in order to increase the stiffness of the aircraft structure, 
a minor modification was made to the finite element model and the 
modulus of the aluminum material was increased without any 
considerations of the other physical characteristics (i.e., specific 
weight, strength allowables, etc.). The modification of the finite 
element model involved only a change of the spar cap properties 
from the aluminum material to a generic graphite-epoxy system 
whose properties were available in the materials data base file. For 
the material modifications, the major modulus of elasticity (En) 
for the Hercules IM7/8552 graphite-epoxy system replaced the 
modulus of elasticity for the aluminum and, in another case, a 
"future improved material" with the modulus of elasticity of steel 
was used. In both cases, the Poisson’s ratio and specific weight of 
aluminum was retained. The pertinent characteristics for each of 
these materials are presented in Table 1. The EZDESIT routine 
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requires that the properties of each of these materials be contained 
in a file called "matprop.prn" which resides in the execution 
directory and contains the data of Table 1 along with other physical 
characteristics such as strength allowables. Since the material 
properties are temperature dependent, a temperature is assigned to 
each element (for the present study, the entire vehicle is assumed to 
be at room temperature, 72 deg. F ) using the PATRAN "DFEG" 
command given on page 22-12 of Reference 2. 

The section properties for each of the elements in the finite 
element model represent the actual structural component (i.e., 
idealization of an integrally stiffened panel as a plate element). The 
lower bound of these section properties is dictated by the minimum 
gage requirements for fabricating and assembling the structural 
components. For the present model, the characteristic section 
properties for the elements used to idealize the structural 
components are defined as follows: 

The rod elements representing the rib and spar caps and the 
longerons of the fuselage are sized by their cross-sectional 
area. 

The beam elements representing the frames of the fuselage are 
sized by the cap or flange area, web height, and web thickness. 

The plate elements representing the skins or surface panels of 
the wing and fuselage are idealizations of honeycomb panels 
and are sized by facesheet thickness and core height. 

The plate elements representing the webs of the ribs and spars 
are idealizations of corrugated panels and are sized by the 
thickness of the material. 

The minimum gage values for the elements used in the present model 
are listed below. 

Rod Area = 0.01 sq. in. 

Beam Cap Area = 0.5 sq. in. 

Web Height = 2.0 in. 

Web Thickness = 0.1 in. 
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Honeycomb Panel 


Facesheet Thickness = 0.01 in. 

Minimum Core Height = 0.1 in. 

Maximum Core Height = 3.0 in. 

Corrugated Panel Web Thickness = 0.01 in. 

To size the elements with EZDESIT, each unique set of finite 
elements (i.e., those elements fabricated in the same way and of the 
same material, such as an aluminum corrugated shear web) must be 
assigned a set of physical property data which specifies the 
necessary design parameters such as the minimum gage sizes shown 
above. These data must be specified during the formulation of the 
model with PATRAN which can be accomplished using the "ADD" 
option of the PATRAN "PROPERTY" command given on page 19-23 of 
Reference 2. The required order of input for these parameters and 
their magnitudes, as used in this study, are excerpted from 
Reference 1 and presented in Appendix C. 

(d) Distribution of Non-S tructural Weight 

The defined weight properties are primarily based on 
interpretations of historical data as it may relate to the HSCT. The 
magnitudes of the non-structural weights and locations of their 
center of gravity is given in Table II of Reference 4 and has been 
reproduced as Table B-2 in Appendix B. The distribution of this non- 
structural weight was derived from an earlier drawing of the 
internal arrangement of the primary components (fuel tanks, 
passenger seating, landing gears, etc.) as reproduced in Figure 3. 
Although the maximum gross weight of Reference 4 differs from the 
flight weights corresponding to the design load conditions given by 
Table A-2 of Appendix A, the ratio of the non-structural weights to 
the maximum gross weight, as given in Reference 4, has been 
retained. Using this ratio to obtain the flight design gross weight 
for a particular load case was expedient, but it did not address the 
issue of weight distribution to assure a trimmed or balanced 
airplane and was probably an unrealistic ratio for most weight 
items except fuel. 

A weighting routine was used in the present study to 
distribute the weight of a particular non-structural item to those 
nodes of the finite element model within the vicinity of the item's 
center of gravity. The affected nodes were visually selected using 
the layout of Figure 3 and/or the center of gravity position given in 
Table II of Reference 4 (Table B-2 of Appendix B). For the selected 
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nodes, the routine calculates and assigns a larger portion of a 
particular item's weight to those nodes nearest its center of gravity 
with a proportional decrease in the amount of weight assigned to 
those nodes further away. The equation used to obtain the weight 
distribution for the non-structural items is as follows: 

Wi- <W/n)[1 - (I, - l av )/(l max - 'min)] 


weight of item assigned to node i, 

total weight of item, 
number of effected nodes, 
distance of node i from CG in X-Y plane, 

(1/n) I lj 
i=1 

the distance of the node furtherest from 
the CG, 

the distance of the node closest to the CG. 

These nodal weights are assigned to the model using the PATRAN 
"DFEG" command for associating temperatures (scalars) with nodes 
(Ref. 2, page 22-12). The EZDESIT routine interprets the 
corresponding data in the PATRAN neutral file as concentrated 
weights. 

(e) "Remeshinq" of Aerod ynamic Pressures 

A computer routine was developed for converting the 
aerodynamic pressure data into a PATRAN neutral file describing the 
nodes and mesh of the aerodynamic model and the magnitude of the 
pressure acting at each node (nodal pressures). PATRAN requires 
that this neutral file be converted from a ASCII text file to a binary 
file. A PATRAN utility program called "reader" (Ref. 2, page 1-23) 
can be used for this conversion. The pressure distribution 
corresponding to a load case superimposed on the mesh of the 
aerodynamic model is shown in Figure 4. 

A PATRAN routine for mapping nodal temperature constraints 
in a global-to-local (or vice versa) application was utilized (Ref. 2, 
page 22-25a) to map the pressure data of the various load cases 
onto the structural mesh of the finite element model. One of this 
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routine’s interactive request for input is the name of the PATRAN 
neutral file which is the name of the binary file discussed above. 
PATRAN recognizes these pressure data as nodal temperature data 
and they are reflected as such in the current PATRAN neutral file 
that is generated. The temperature data are converted into pressure 
data during the EZDESIT execution which is described below. 

(f) EZDESIT Requirements 

Since EZDESIT accepts only element temperature data as input, 
the nodal temperatures contained in the PATRAN neutral file, 
discussed above, were converted to element temperatures using an 
interactive program developed for this purpose. Options contained in 
the program allow for assigning element temperatures to any or all 
of the named components discussed in Section 3.0 (a). This program 
then creates a new neutral file which must be called "patran.out" due 
to EZDESIT input specifications. 

4.0 EZDESIT Execution 


(a) Generation of the COMET Runstream 

EZDESIT is a menu driven program with interactive queries 
which culminates into a COMET runstream. The first user response, 
prior to the appearance of the menu, is for the number of bars, 
beams, triangular plates, and quadrilateral plates (in this order) 
which is 

577 229 32 476 

for the present finite element model. The program's main menu then 
appears on the screen and is reproduced below, for clarity, since it 
is similar, but not identical, to that menu documented in Reference 
2 . 
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1 - Create a COMET runstream through XQT INV 

2 - Not available 

3 - Create a partial COMET runstream 

for applying/combining loads and constraints 

4 - Update COMET element stiffness matrices 

5 - Translate COMET results to PATRAN results files 

6 - Run element sizing code 

7 - Create worst case dimension file 

8 - Review calculated component weights 

9 - Exit to system 

The initial selection of menu item 1 translates the nodal 
locations, element connectivity, constraint cases, and nodal weights 
from the PATRAN neutral file (patran.out) into a COMET runstream. 
Also, dummy data are inserted into the runstream for the material 
and section properties and nonstructural weights. The so-called 
"nonstructural weights" represent the weight of the individual 
structural finite elements. Next, menu item 6 is selected to create 
the initial unit weight file, called "uwtmin.out", and dimension file, 
called "Idstmin.dim", for the elements which are based on the 
specified material(s) and minimum gage sizing contained in the 
PATRAN neutral file. The material data are contained in a file called 
"matprop.prn". Menu item 4 is then selected to create initial 
element stiffnesses based on the dimension files generated by 
execution of menu item 6 and the physical characteristics of 
corresponding materials which are used to update the dummy data 
created under menu item 1. These element stiffnesses are generated 
by a program called STIMAT which has been incorporated as a 
subroutine in the EZDESIT procedure. Finally, selection of menu item 
3 converts PATRAN applied load information to a COMET applied load 
format, sets up inertial load cases described by the analyst, and 
combines several load cases to define a design load condition. This 
information and data are then assembled into a portion of the COMET 
runstream that can include processors for generation of stress 
resultants and nodal displacements files. Selection of menu item 9 
returns the user to the operating system. 

At this point, it is necessary to edit the runstream file 
created by menu item 3 to accomplish the following: 

Remove the system commands contained in the first four lines 

of the COMET runstream. 
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The following statement is included immediately after entry 
into the E processor since most quadrilateral plate elements 
used for modeling the the wing and fuselage cover panels 
exceed the default element warping (i.e., not all four nodes in a 
plane) limits of the COMET program. 

T= 1.-20, .05, 1.-5, 1.-1, 20., 1.-4, 1.-4, 1.-4 

Although excessive warping does not produce a fatal error, a 
notice of failure is written into the log file at each occurrence 
and the revision of the fourth term from its default value of 
I.OxlO' 5 to the indicated value of I.OxlO' 1 will substantially 

reduce this output. 

Insert "eldata: pres exx sid 1" for the present study at the 
appropriate locations where PATRAN element temperature data 
were used to represent element pressure data. (The "xx" of exx 
is the element number such as e43 or e33 and "sid" is the load 

case number.) 

Due to an EZDESIT formatting error, the "summing" statement 
for combining the aerodynamic pressures and the inertial loads 

requires editing. 

Move to the end of the runstream file and insert "[XQT EXIT". 

This edited version of the loads and solution portion of the COMET 
runstream is merged at the end of the runstream created by menu 
item 1 for a complete runstream that can be submitted for batch 
execution of the element sizing process. An abbreviated listing of a 
runstream for the all-titanium finite element model is presented in 
Appendix D. As will be indicated in the next section, this runstream 

(file) is called "tb.dat". 


(b) Ratch Sub m ittal of the COMET Runstream 

The EZDESIT program creates a COMET runstream, but the 
iterative procedure for resizing the structure is accomplished using 
a program called EZBATCH. The EZBATCH program accesses the 
COMET system and submits the COMET runstream (tb.dat) by calling a 
command (script) file called "tb.scr". (If the name of the runstream 
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is other than "tb.dat", "tb.scr" must be modified to reflect this new 
name.) A listing of "tb.scr" follows: 


#/!bin/csh -f 
source /csm/login 
testbed ctb.dat 

Another file called "ezbatch.dat" must be created. This file defines 
the load sets, constraint sets, associated thermal load sets, and 
associated element temperatures that are used for the iterative 
resizing process. A description of the contents of this file is given 
in Reference 2 and a listing of the file as it relates to the present 
study follows: 


577 229 32 476 
5 

51 52 53 61 62 
11111 
0 0 0 0 0 
0 0 0 0 0 
0 0 0 0 0 
11111 
5 

Next, a copy of the file called "Idstmin.dim" (created during 
execution of menu item 6) is copied to a file named 
"Idstcomb.dimjOO" which is used for creating a worst-case element 
dimension file. Finally, the program "EZBATCH" is submitted as a 
batch job to initiate the element sizing process. A typical UNIX 
command for submitting the job is 

ezbatch >& ezbatch.out & 

which assumes the program execution file resides in the current 
working directory. The file "ezbatch.out" contains iteration 
information and is created as the EZBATCH execution proceeds. 



5.0 Design Cycle 

(a) load Iterations for B alanced Airplane 

For the initial structural resizing, the airplane was not 
balanced for free flight which usually results in large boundary 
reactions at the fictitious boundary constraints. This imbalance 
was unavoidable due to the unavailabilty of the true structural 
weight prior to structural sizing and, consequently, the magnitude of 
the resultant lift component required to balance and trim the 
airplane. Once the structure was resized for all five design load 
conditions, an iterative process was carried out to balance the 
aerodynamic pressure loads with the inertial forces (structural and 
non-structural weights). A control force was applied at the aft end 
of the airplane to counteract the moment caused by the difference in 
the center of gravity and the center of pressure for each of the five 
design load cases. Usually, two iterations were required to obtain a 
nearly balanced airplane (which minimized the boundary reactions) 
after the resultants of the aerodynamic pressure loads and the 
inertial forces were balanced. The lift component of the 
aerodynamic load was increased or decreased to compensate for the 
vertical balance load. A detailed description of this balancing 
process is given in Appendix E for a typical design load condition. 

(b) Resizing Iterations 

In the course of obtaining a balanced airplane for the initially 
sized structural elements, it was necessary to modify the applied 
loads and weights. For the present study, a second resizing of the 
structural elements was performed only on the all-titanium finite 
element model. This resizing involved another execution of the 
EZDESIT program to incorporate the new applied forces and to create 
a new COMET runstream. However, the experienced user can avoid a 
complete re-execution of EZDESIT by carefully editing the COMET 
runstream and relocating or renaming the output files from the 
initial EZBATCH execution. 

6.0 Results 

The initial structural resizing results for the all-aluminum 
finite element model of the candidate Mach 3.0 HSCT yielded large 
wing tip displacements. These large displacements are primarily 
due to the small depth of the external shape of the wing and 
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secondarily due to the small number of spars and their arrangement 
within the wing. Variations of the external shape or the structural 
arrangement to obtain increased stiffness was not done for this 
initial investigation because of the excessive time required to 
generate a new finite element model. However, variations of the 
structural stiffness were obtained by using different materials. The 
various non-metallic materials considered were based entirely on 
the major modulus of elasticity of the composite laminates and the 
resizing algorithm was based on the failure criteria and 
manufacturing requirements (i.e., minimum gage) of metallic 
materials and their associated allowable design strengths. 

(a) Weight 

The structural weights of the finite element models for the 
various materials considered are shown in Table 2. These structural 
weights do not include a significant portion of the structural weight 
due to fasteners, splices, shims, etc. A second resizing iteration of 
the all-titanium model resulted in a structural weight increase of 
nearly 13%. This weight increase is probably due to the increase in 
lift on the wing that was required to obtain a balanced airplane and 
resulted in higher internal loads, particularly at the wing-fuselage 
intersection. No aerodynamic forces were applied to the fuselage 
and the aerodynamic forces were applied to the wing's lower 
surface. A third iteration resulted in a negligible change of this 
weight. 

(b) Stiffness and Deflection 

A comparison of the maximum wing tip deflection (at the 
trailing edge of the tip rib) corresponding to the various design load 
cases and the finite element models composed of various materials 
is presented in Table 3. The resizing algorithm did not include a 
deflection (or stiffness) constraint which resulted in the large wing 
tip deflections for this primarily strength designed structure with 
minimum gage constraints. Of course, increasing the elastic moduli 
of the materials will produce smaller deflections. But stiffening 
the spar caps using a graphite-epoxy material for an otherwise all- 
aluminum model resulted in little or no change in the deflections. 

The ineffectiveness of the stiffened spar caps is due to the fact that 
the overall stiffness effects (i.e., EA, El, and GJ) of the wing skins 
are so much greater than the stiffness contributed by the spar caps. 
There was no variation of the minimum gages for the present study. 
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The deflections resulting from another resizing iteration of 
the all-titanium model with the new balance loads indicated a 
decrease in the tip deflection of approximately 30%. This decrease 
in tip deflection was probably due to the decreased wing loading 
(caused by compensating for the balance load) and the relieving 
effect of the increased structural weight. Another resizing 
iteration for this model showed little, if any, change in the wing tip 
displacement. 

The quadrilateral elements representing the cover panels of 
the wing and fuselage are excessively warped with respect to the 
element specifications of COMET. The effect of this warping on the 
shape of the deformed structure is not known; however, results 
presented in Reference 7 indicated small but significant effects (4% 
to 7%) on the displacement of a simple, cantilevered, box-beam 
structure. 

7.0 Concluding Remarks 

A structural sizing study for a single configuration of a 
candidate Mach 3.0 High-Speed Civil Transport has been conducted. 

An integrated analysis system incorporating routines to generate 
structural arrangements (and corresponding finite element models) 
along with the routines used herein provide a design tool that is both 
timely and cost effective for performing such analyses. In this 
regard, a detailed description of the procedures used and the data 
handling tasks performed during the execution of this study have 
been included in this report to facilitate the implementation of a 
structural sizing procedure into such a system. 

The wing tip displacements are excessive for all materials 
considered for the finite element model used in the study. Major 
increases in stiffness will only be obtained by increasing the depth 
and/or rearranging the internal structure (ribs and spars) of the 
wing. Significant increases in stiffness can be obtained by using 
structural components made of composite materials; however, the 
present version of EZDESIT is limited to the resizing of structural 
components made of isotropic materials. Furthermore, the increased 
strength of the composite materials would result in a lighter 
vehicle; hence, reduced loads. Finally, refinement of the 
displacements can be obtained by evaluating the effect of the 
warped plate elements and mesh refinements. 
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Table 1 - Physical characteristics of the materials 
used in the present study. 


Material 

Modulus 

of 

Elasticity, 

msi 

Poisson's 

Ratio 

Specific 

Weight, 

3 

Ibs/in. 

2024-T81 

Aluminum 

10.0 

0.33 

0.100 

Ti-6-4 

Titanium 

16.0 

0.31 

0.160 

IM7/8552 * 
Graphite-Epoxy 

24.5 

0.33 

0.100 

Generic 

Graphite-Epoxy 

19.3 

0.30 

0.058 

Future 

Improved 

Material 

30.0 

0.33 

0.100 


* Product of Hercules Materials & Systems Company 
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Table 2 - Comparison of the structural weight for 
finite element models composed of various materials. 


Material 

Structural Weight, lbs 

Design Iteration 1 

Titanium 

34,256.8 

Aluminum 

35,871.8 

IM7/8552 

32,272.2 

Future Improved 

31,680.6 

Material 


Aluminum with 

35,812.0 

Graphite-Epoxy 


Spar Caps 


Design Iteration 2 

Titanium 

38,657.4 

Design Iteration 3 

Titanium 

38,684.4 















Table 3 - Comparison of the maximum wing tip displacement for 
the design load cases with finite element models composed of 

various materials. 


cjnjfp Figment 

Wing Tip Displacements, inches 

I Iflllt? t-lwl 1 Ivl 1 1 

Model 

Case 1 

Case 2 

Case 3 

Case 4 

Case 5 


Design Iteration 1 

Titanium 

220.0 

251.0 

144.0 

640.0 

667.0 

Aluminum 

179.0 

199.0 

116.0 

520.0 

537.0 

IM7/8552 

76.7 

84.7 

49.4 

223.0 

229.0 

Future Improved 

63.5 

70.1 

40.9 

184.0 

190.0 

Material 






Aluminum with 
Graphite-Epoxy 

180.0 

199.0 

116.0 

521.0 

537.0 

Spar Caps 







Design Iteration 2 

Titanium 

158.0 

177.0 

102.0 

458.0 

472.0 


Design Iteration 3 

Titanium 

158.0 

176.0 

102.0 

457.0 

471.0 
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Figure 1. General arrange 



20 


Figure 2. Structural arrangement and finite element model. 



21 


Figure 3. Internal arrangement of primary non-structural weight items. 
The maximum fuel weight (lbs.) of the various tanks is indicated. 






Appendix A 


Weight Summary and Design Flight Loads for a Candidate 
Mach 3.0 High-Speed Civil Transport 

The baseline flight load cases used for determining the 
aerodynamic and inertia loadings for the present study are given in 
Table A-1. These flight conditions were based on trajectory 
optimization studies for a version of the HSCT whose weight 
breakdown is presented in Table A-2. 
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Table A-1 - Baseline flight load cases for a candidate 

Mach 3.0 HSCT. 


Load 

Case 

Mach 

Number 

Altitude, 

ft. 

Load 

Factor, 

9's 

Flight 

Weight, 

lbs 

C 

L 

Comments 

1 

3.0 

72,700 

1.0 

494,000 

0.077 

Mid cruise 

2 

1.2 

21,300 

1.0 

656,000 

0.058 

Transonic climb 

3 

■ 

43,000 

1.0 

333,000 

0.141 

Reserve cruise 

4 

3.0 

59,000 

2.5 

623,000 

0.127 

Structural limit at 
cruise Mach number 

5 

0.6 

10,000 

2.5 

664,000 

0.363 

Structural limit at 
low speed 
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Table A-2 - Weight summary for a version of a candidate 

Mach 3.0 HSCT. 


ITEM WE I GHT, l b s 


Wing 

Horizontal Tail 
Vertical Tail 
Canard 
Fuselage 
Landing Gear 
Nacelles 

Total Structure 

Engines 

Thrust Reverser 
Propeller 
Gear Box 

Miscellaneous Systems 

Plumbing 

Tanks 

Insulation 

Total Propulsion 

Surface Controls 

Auxiliary Power 

Instruments 

Hydraulics 

Electrical 

Avionics 

Furnishings and Equipment 
Air Conditioning 
Anti-icing 

Total Systems and Equipment 


76974 

0 

554 

0 

43906 

22615 

18440 

162489 

23096 

0 

0 

0 

1508 

5694 

0 

0 

30298 

9240 

1210 

1410 

5526 

4170 

1732 

20302 

7827 

274 

51691 
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Table A-2 - Con't. 


ITEM WEIGHT, lbs 


EMPTY WEIGHT 


244478 

Flight Crew 

450 


Cabin Crew 

1130 


Unusable Fuel 

3636 


Engine Oil 

545 


Passenger Service 

3560 


Cargo Containers 

7162 

16483 

Total Operating Weight 


260961 

OPERATING EMPTY WEIGHT 

Passengers (250) 

41250 


Baggage 

11000 


Cargo 

2565 


Total Payload 


54815 

Zero Fuel Weight 


315776 

Mission Fuel Weight 


359924 

MAXIMUM GROSS WEIGHT 


675700 
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Appendix B 


External Shape and Mass Properties for a Baseline Mach 3.0 

HSCT Concept 

The finite element model used in the present study was 
generated within the external shape of a baseline HSCT concept 
defined in Reference 4 and is reproduced herein as Table B-1. The 
distribution of the non-structural mass for this finite element 
model was based on empirically estimated weight and balance data 
for this baseline concept as given in Table B-2 
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Table B-1 - Numerical model for the external shape of a 
baseline Mach 3.0 HSCT concept. 


1 

12183 * 

0.0 

23.0 

0.0 

1.163 

7.270 


2 A 8 T 3 I 
1 0 1 
129.70 
.30 

30.0 

0.0 

2.0 

3.0 


FUSELAGE AS WING 
0 0 2 14 20 

160.0 


ANL2/SU0AERF MODEL including^pqds 


.73 

1.23 2 

33.0 

40.0 S ' 

0.0 

209. 147 

-.70 

207.984 

- 1.70 

201.877 

- 2.60 

192.789 

- 3.37 

183.814 

- 4.80 

133 . 0 S 2 

- 6.69 

1 23. 909 

- 8.28 

91.110 

- 9.00 

33.309 

- 0.00 

47 . 2 B 1 

- 6 . 47 

36.723 

- 3.60 

25 . 130 

- 3.90 

21.320 

- 6.00 

8.400 


30 

>.0 


3.00 

60.0 


7.30 

70.0 


10. 0 

80.0 


13.0 

90.0 


20.0 

100.0 


16.338 7.3 
23.333 9.333 
36.096 13.00 
83.238 20.00 

121 . 23827.00 
160 . 33334.202 

171 . 33340.00 

188 . 66733.00 

210 . 49368.00 

219 . 63872.00 
238.40 73.30 - 

0.000 - 0.069 - 0.102 - 0.176 
- 3.700 - 4.230 - 4.730 - 3.196 
0. 000 - 0.069 - 0.102 - 0.176 
- 3.927 - 4.480 - 4.939 - 3.366 
0.000 - 0.028 - 0.044 - 0.079 
- 3.683 - 4.420 - 3.009 - 3.407 
0.000 0.027 0.037 0.032 

- 3.138 - 3.908 - 4.623 - 3.273 
0.000 0.031 0.048 0.073 

- 2.209 - 2.917 - 3.334 - 4.132 
0.000 0.064 0.096 0.133 

- 0,467 - 0.943 - 1.442 - 1.939 
0.000 0.063 0.097 0.137 

- 0.206 - 0.340 
0 . 07 B 0.128 


SCXCS 

XAF1 

XAF2 

WORG 

WORG 

WORG 

WORG 

WORG 

WORG 

WORG 

WORG 

WORG 

WORG 

WORG 

WORG 

WORG 

WORG 


1 

1A 


4 

3 

6 

7 

8 
9 

10 

11 

12 

13 


0.338 

0.000 

0.887 

0.000 

0.049 

0.000 

0.077 


0 . 103 
0.033 
0.831 
0.023 
0 . 682 
0.002 
0.089 


0.730 
0.034 
0.704 
0.003 
0 . 103 


0.606 
0.036 
0.726 
0.003 
0. 114 


0.000 - 0.006 
- 0.233 - 0.293 
0.000 - 0.002 
- 0.079 - 0.093 
0.000 0.000 
- 0.004 - 0.006 
0.000 0.000 
0.008 0.009 


-0.009 -0.013 
-0.333 -0.372 
-0.003 -0.004 
-0.107 -0.120 
0.000 0.000 
-0.008 -0.011 
0.000 0.000 
0.011 0.013 


- 0.360 

- 3.933 

— 0 • s>68 

- 6.073 
- O . 190 
- 6.241 
0.032 
- 6.292 
0 . 160 
- 3. 163 
0.286 
- 2.983 
0.293 
- 1.260 
0.247 
0.292 
0 . 112 
0.764 
0.011 
0 . 134 
- 0.029 
- 0.437 

-o.ooa 

- 0.147 

0.000 

- 0.020 

0.001 

0.016 


-0.811 
-6.383 
-0.81 1 
-6.680 
-0.480 
-6.800 
-0. 130 
-7. 133 
0.217 
-3.968 
0.484 
-3.946 
0.310 
- 2.001 
0.432 
-0.094 
0.218 
0.786 
0.021 
0 . 144 
-0.037 
-0.487 
-0.017 
-0. 168 
0 . 000 
>0.027 
0.002 
0.019 


-1.261 
-7. 160 
- 1.262 
-7.226 
-0.832 
-7.234 
-0. 147 
-7.773 
0. 037 
-6.603 
0.343 
-4.803 
0.638 
-2.722 

0.614 
-0.307 
0.320 
0.806 
0.030 
0. 130 
-0.084 
-0.321 
-0.023 
- 0. 181 
0.000 
-0.036 
0.002 
0.022 


-1.699 
-7.704 
-1.718 
-7.7*1 
- 1.211 
-7.671 
-0.767 
-8.239 
- 0.202 
-7.082 
0.304 
-3.323 
0.700 
-3.393 
0.723 
-0.923 
0.418 
0.324 
0.03? 
0 . 1*9 
-o.lll 
-0.337 
-0.032 
-0. 189 
0.000 
-0.0*3 
0.003 
0.023 


-2. 478 
-9.219 
-2.360 
— 0 • 2<j0 
-2.039 
-8.079 
-1.373 
-9.333 
-0.798 
-7.436 
0.300 
-6.087 
0.690 
-3.990 
0.841 
-1.326 
0.330 
0.827 
0.034 
0. 131 
-O. 162 
-0.344 

-o.o*a 
-0. 190 
-0. 001 
-0. 044 
0.005 
0.023 


-3. 117 

ZORD t-1 

-8.716 

ZORD 1-2 

-3.288 

Z0RD1A-1 

-8.731 

Z0RD1A-2 

-2.393 

ZORD 2-1 

-8.463 

ZORD 2-2 

-2.392 

ZORD 3-1 

-8.731 

ZORD 3-2 


-1.487 
-7.73? 
-0.049 
-6.460 
0.370 
-4. 473 
0.380 
-1. 697 

0.6H4 

0.330 
0.066 
0 . 140 
-0 . 209 
-0.348 
-0.063 
-0. 187 
- 0. 002 
- 0 . 0*0 
0.006 
0.031 


ZORD 4-2 
ZORD 5-1 
ZORD 3-2 
ZORD 6-1 
ZORD 6-2 
ZORD 7-1 
ZORD 7-2 
ZORD 3-1 
ZORD 3-: 
ZORD 9-1 
ZORD 9-2 
ZORD 10- 1 
ZORD 10-2 
ZORD 11-1 
ZORD 11-2 
ZORD 1 
Z0RD1 
ZORD 1 
ZORD 1 
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(I r i n t> 


Table B-1 - Con’t. 


0.0 

. 10 

. 14 

^ T 

2.23 

2.41 

2.36 

2.68 

0.0 

. 08 

. 14 

.20 

1.96 

2. 17 

2 • 

2. 49 

0.0 

.08 

. 12 

.21 

1.08 

1.31 

1 . SB 

1.93 

0.0 

.061 

. 100 

. 102 

.886 

.953 

.990 

.996 

0.0 

. 1333 

. 1399 

.2037 

.B4JZ 

.8987 

.9330 

.9311 

0.0 

. 1424 

. 1709 

.2177 

.9034 

.9603 

.9970 

1.0163 

0.0 

. 1373 

. 1887 

.2404 

.9999 

1 . 0603 

1. 1010 

1 . 1223 

0.0 

.2013 

.2429 

.3087 

1 . 2088 

1.2801 

1.3280 

1 . 3323 

0.0 

.0313 

.0769 

. 1289 

L 9 * 

1.8341 

1.9019 

1 . 9330 

0.0 

.0332 

.0793 

. 1330 

1.7880 

1 . 0920 

1 . 9620 

1.9970 

0.0 

. 0307 

.0737 

. 1269 

1.7032 

1.8044 

1.8712 

1.9043 

0.0 

.0463 

.0694 

. 1 174 

1.5643 

1 . 6333 

1.7166 

1.7474 

0.0 

. 0328 

.0709 

. 1323 

1 . 738 Q 

1.8600 

1.9272 

1.9616 

C.O 

.0328 

.0709 

. 1323 

1 . 7589 

1 . 8600 

1.9272 

1 . 9616 

180.0 

9.7 

-14. 70 


0.0 

5.0 

10.0 

15.0 

3. 193 

3.307 

• 336 

3.670 

180.0 

16. 7 

-14.27 


0.0 

3.0 

10.0 

13.0 

3. 195 

3.387 

3.336 

3.670 

180.0 

23.7 

-13.29 


0.0 

5.0 

10.0 

15.0 

3. 193 

3.387 

3.336 

3.670 

171.75 

11.4 

-12.00 


0.0 

3.25 

5.73 

8.23 

0.0 

.626 

.937 

1 . 12 B 

171.73 

21.8 

-11.0 


0.0 

3.25 

3.73 

9.23 

0.0 

.626 

.937 

1.128 


. 44 

.7* 

1 . 07 

1.30 

2.93 

2.97 

3. 43 

3.70 

.37 

.62 

.70 

. 92 

2.66 

2.93 

3. 17 

3 . 29 

.30 

.41 

.30 

.59 

*** *T— ’ 

■ w w 

2 . 29 

1.97 

1.68 

.287 

.405 

.480 

.362 

.994 

.962 

.838 

.586 

.2970 

.4131 

. 3017 

.3741 

.9330 

.8340 

.6731 

.4636 

.3174 

.4433 

.3631 

.6134 

.9970 

.8912 

.7192 

.4973 

.3303 

.4898 

.3921 

.6774 

1. 1010 

.9841 

.7943 

.3494 

. 4324 

.5952 

.7139 

.8197 

1.3244 

1. 1826 

.9333 

.6581 

.2327 

.4939 

.7314 

.9337 

1.0942 

1.6896 

1.3591 

.9374 

.2607 

.3093 

.7345 

.9839 

1 . 9340 

1.7430 

1 . 4020 

.9672 

.2486 

. 4839 

.7196 

.9402 

1 . 8633 

1 . 6623 

1.3371 

.9222 

.2281 

. 4438 

• 6602 

.8627 

1.7099 

1.3231 

1.2268 

.8461 

.2376 

.4997 

.7411 

.9690 

1.9110 

1.6936 

1.3356 

.9312 

.2376 

.4997 

.7411 

. 9698 

1.9110 

1 . 6936 

1 . 3336 

.9312 

20.0 

23.0 

27.3 

30.0 

3.810 

3.937 

3.999 

3.9ei 

20.0 

23.0 

27.5 

30.0 

3.810 

3.937 

3. 989 

3.901 

20.0 

23.0 

27.3 

30.0 

3.810 

3.937 

3.989 

3.901 

10.75 

13.23 

13.73 

18.23 

1.223 

1.194 

1.095 

.982 

10.73 

13.25 

13.73 

18.23 

1 . 223 . 

1. 194 

1.093 

. 982 


1.72 

2.03 

WORD 1-1 

3.35 

2.69 

WORD 1-2 

1.34 

1.70 

WORT 1 A 1 

2. 92 


WORD 1A2. 

. 70 

.86 

WORD 2-1 

1.41 

1.04 

WORD 2-2 

.685 

.793 

WORD 3-1 

.283 

0.0 

WORD 3-2 

.6893 

.7708 

WORD 4-1 

.2352 

0 . 0 

WORD 4-2 

. 7365 

.8322 

WORD 5-1 

.2513 

0 . 0 

WORD 5-2 

.8134 

.9190 

WORD 6-1 

.2775 

0.0 

WORD 6-2 

. 9049 

1.1113 

WORD 7-1 

.3331 

0 . 0 

WORD 7-2 

1.3402 

1.3946 

WORD 8-1 

.4750 

0 . 0 

WORD 8-2 

1.3908 

1.6430 

WORD 9-1 

. 4900 

0.0 

WORD °-2 

1 . 3264 

1.5688 

WORD 10-1 

. 4673 

0.0 

WORD 10-2 

1.2169 

1.4394 

WORD 11-1 

.4288 

0.0 

WORD 1 1-2 

1 . 3683 

1.6196 

WORD 12-1 

. 4710 

0.0 

WORD 12-2 

1.3605 

1 .6196 

WORD 13- 1 

.4710 

0.0 

WORD 1 3-2 
PODORG 1 

33.0 

38. 47 

XPOD 1 

3.093 

3.004 

PCDR 1 
PODORG 2 

33.0 

38. 47 

XPOD 2 

3.093 

3.804 

PQDR 2 
POGORG 3 

33.0 

30.47 

XPOD 3 

3.893 

3.004 

PODR 3 
PODORG 4 

23. 23 

20. 2S k 

XPOD 4 

• 67 

0.0 

PODR 4 
PODORG 3 

23.23 

29.23 

XPOD 3 

.67 

0.0 

PODR 3 
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Table B-2 - Estimated weights and balance for a 
baseline Mach 3.0 HSCT concept. 


Item 


Wing 

Horizontal ■ . 

Vortical ... 

Fin 

Canard .... 
Fuselage 
Liinding gear 
N;icelle 

Structure lotal 


Weight, lb 

■ 1 

FS c.g.. in 

WL c.g.. 

j 56258 

2050.0 

85.0 

i 

0 


; 

! t ioi 

3462.1 

297.3 

! ° 



0 



43906 

1920.0 

103.0 

! 32326 

1850.4 

70.0 

! 19 790 

2330.0 

13.0 

j 153381 

2010.9 

j 79.2 


Engines 

Thrust reverwer 
Miscellaneous systems 
Fuel system . . 
Propulsion total 


50665 

0 

2341 
8 026 
61031 

2435.0 

1481.5 

2083.9 

2352.2 

13.0 

45.7 

89.5 

24.3 

7 399 

2193.6 

122.2 

1578 

2600.0 

50.0 

2 141 

1172.4 

103.0 

5 449 

1958.G 

103.0 

3888 

1293.5 

103.0 

1 449 

860.0 

136.0 

20 388 

1842.0 


6 741 

1890.0 

■ 

m 

1948.0 

654 

49 246 

1837.7 

61.2 


Surface controls 
Auxiliary power unit 
Instruments 
Hydraulics 

Electrical 

Avionics 

Furnishings and equipment 
Air conditioning 
Anti-icing 


System ana equipment 

— ! 

Weight empty 

263 658 


63.1 

i 

| Flight crew and baggage (2) 

450 

528.0 

144.0 j 

98.0 \ 

Cabin crew and baggage (7) . 

1 130 

1842.0 

i 1 

* { tiusuable fuel j 

3 269 

2083.9 

89.5 j 

k ion 

9M 

2435.0 

13.0 ! 

Engine oil : 

Passenger service j 

3 560 
0 

1842.0 

98.0 

- — 

Operating weight 

272981 

2056.3 

64.0 



Passengers (250) 

41250 

1842.0 

! 98.0 1 

1 1 1 000 

1300.0 


Passenger baggage 

\ 2 565 

1920.0 

103.0 

j Miscellaneous items 

f’.irffn 

i 0 


t 

Zero-fuel weight 

; 327 796 

2002.9 

69.0 

' 



J Mission fuel 

j 385900 

1649.0 

106.1 

| (bos* weight 

713696 

1919.7 

j 89.0 
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Appendix C 


Description of Physical Property Data Requirements for an 

EZDESIT Execution 

Tables C-1 through C-4 list the existing property set 
definitions for types of element constructions that can currently be 
used with EZDESIT. The format and description of the data is taken 
verbatim from Reference 2 and the value used for each of these 
items in the present study is listed. 
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Table C-1 - Physical property definition for an isotropic 

material, honeycomb plate. 


FIELD CURRENT VALUE DESCRIPTION 


1 

10 a 

Material number, refer to the 
"matprop.prn" file 

2 

30 

Span used to check panel 
buckling, in. 

3 

9 

9 = isotropic honeycomb 

4 

1.0 

Factor times input loads to 
obtain ultimate load (typically 
1 .0 with EAL runs done at 
ultimate load) 

5 

0.8 

Panel buckling load knockdown 
factor (typically .8) 

6 

1.0 

Non optimum factor (1.0 = 
none) 

7 

0.01 

Minimum gage of facesheet, in. 

8 

1.0 

Young’s modulus reduction 
factor (1.0 = full value in 
"matprop.prn" file) 

9 

1.0 

Allowable strength reduction 
factor (1.0 = full value in 
"matprop.prn" file) 

10 

0 

"Dummy" b 

1 1 

0.1 

Minimum core height, in. 

12 

3.0 

Maximum core height, in. 

13 

0.1 

Core material density, lbs/in. 3 


(ex: aluminum = .1, routine 
hcbndth assumes 2% core 
factor) 


14 

0 

"Dummy" 

15 

0 

"Dummy" 

16 

0 

"Dummy" 

17 

0 

"Dummy" 

18 

0 

"Dummy" 

19 

0 

"Dummy" 

20 

0 

"Dummy" 


a This "material number" is for titanium. 

b The "Dummy" parameters are only placeholders for possible code expansion. 
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Table C-1 ■ Con’t. 


HELD 

CURRENT VALUE 

DESCRIPTION 

21 

0 

"Dummy" 

22 

0 

"Dummy" 

23 

0 

"Dummy" 

24 

0 

"Dummy" 

25 

0 

"Dummy" 

26 

0.67 

Factor times ultimate load to 
obtain limit load level (ex: = 

.67 if using a 1.5 safety factor) 

27 

0 

"Dummy" 

28 

0 

"Dummy" 

29 

0 

"Dummy" 

30 

0 

"Dummy" 

31 

0 

"Dummy" 

32 

0 

"Dummy" 

33 

0 

"Dummy" 

34 

0 

"Dummy" 

35 

0 

"Dummy" 
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Table C-2 - Physical property definition of an isotropic 

material, corrugated web. 


FIELD 


P.l IRRFNT VALUE DESCRIPTION 


1 

2 

3 

4 

5 

6 


7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 


19 

20 
21 
22 
23 


6 a 

30 

10 

1.0 

0.8 

2 


0.01 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0.67 


0 

0 

0 

0 

1.0 


Material number, refer to 
the "matprop.prn" file 
Panel buckling span, in. 

10 = corrugated web 
Factor times input loads 
to obtain ultimate load 

Panel buckling load 
knockdown factor typically .8 
empirical/theoretical ratio) 

= 1 corrugated along element x 
direction 

= 2 corrugated along element y 
direction 

Minimum material gage, in. 

"Dummy" b 

"Dummy" 

"Dummy" 

"Dummy" 

"Dummy" 

"Dummy" 

"Dummy" 

"Dummy" 

"Dummy" 

"Dummy" 

Factor times ultimate load to 
obtain limit load level (ex: *= 

.67 if using 1.5 safety factor) 
"Dummy" 

"Dummy" 

"Dummy" 

"Dummy" 

Modulus reduction factor (1.0 * 
full value from "matprop.prn" 
file) 


a This "material number" is for titanium. 

b jhQ "Dummy" parameters are only placeholders for possible code expansion. 
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Table C-2 - Con’t. 


FIELD 

CURRENT VALUE 

DESCRIPTION 

24 

1.0 

Allowable strength reduction 



factor (1.0 = full value from 



"matprop.prn" file) 

25 

1.0 

Non optimum factor (1.0 = 


none) 
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Table C-3 - Physical property definition for a bar. 


FIELD 

Cl IRRFNT VALUE 

DESCRIPTION 

1 

6 a 

Material number, refer to 
"matprop.prn" file 

2 

1.0 

Non optimum factor 

3 

13 

13 - bar 

4 

1.0 

Factor times input loads to 
obtain ultimate load 

5 

0.01 

Minimum area, in. 2 

6 

1.0 

Modulus reduction factor (1.0 = 
full value from "matprop.prn" 
file) 

7 

1.0 

Allowable strength reduction 
factor (1.0 = full value from 
"matprop.prn" file) 

8 

0 

"Dummy" b 

9 

0 

"Dummy" 

1 0 

0 

"Dummy" 

1 1 

0.67 

Factor time ultimate load to 
obtain limit load level (ex: = . 
67 if using a 1.5 safety factor) 

12 

0 

"Dummy" 

1 3 

0 

"Dummy" 

14 

0 

"Dummy" 

15 

0 

"Dummy' 


a This "material number" is for titanium. 

b The "Dummy" parameters are only placeholders for possible code expansion. 
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Table C-4 - Physical property definition for a planar beam. 


FIELD CURRENT VALUE DESCRIPTION 


1 

2 

3 

4 

5 

6 

7 

8 


9 


10 

11 

12 

13 


14 

15 

16 
17 


6 a 


1.0 

14 

1.0 

0.5 

2.0 

0.1 

1.0 


1.0 


0 

0 

0 

0.67 


0 

0 

0 

O' 


Material number, refer to 
the "matprop.prn" file 

Non optimum factor (1. = 
none) 

14 = planar beam 
Factor times input loads 
to obtain ultimate load 
Minimum cap area, in. 2 
Minimum web height, in. 

Minimum web gage, in. 

Modulus reduction factor 
(1.0 « full value from the 
"matprop.prn" file) 

Allowable strength 
reduction factor (1.0 = full 
value from the "matprop.prn" 
file) 

"Dummy" b 

"Dummy" 

"Dummy" 

Factor time ultimate load to 
obtain limit load level (ex: = 

.67. if using a 1.5 safety factor) 
"Dummy" 

"Dummy" 

"Dummy" 

"Dummy" 


a This "material number" is for titanium. 

b The "Dummy" parameters are only placeholders for possible code expansion. 
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Appendix D 

COMET Runstream for an All-Titanium HSCT 
Finite Element Model 

An abbreviated listing of a COMET runstream, called "tb_dat" 
that was submitted for structural resizing of a candidate Mach 3.0 
all-titanium HSCT finite element model is presented in this 
Appendix. 


♦d»f/i ns_co*»on«200000 

•open 1 ti.101 /new 
•set echo»off 
♦add genutil.prc 
CXOT TAB 
START 398 

TITLE ’ALL TITANIUM HSCT MODEL 7/30/91 4:30 
HATC 

1 15945143. 0.31 0.000 0.481710E-05 0.481710E-05 0. 

2 15945143. 0.31 0.000 0.481710E-05 0.481710E-05 0. 

NSH 

1 0.00167 

2 0.00417 

3 0.00333 


• 

135 

136 

137 
JL0C 


1 

2 

3 


396 

397 

398 


0.00972 

0.01111 

0.02917 


30.00000 

30.00000 

30.00000 


3480.00000 

3480.00000 

3480.00000 


0.00000 

11.33280 

0.00000 


10. 18640 
3.09320 
0.00000 


-7.41960 

-7.19358 

0.82680 


-13.54595 

-10.82704 

- 10.02000 


-1REF 

F0RMAT=2 
1 1 


0.00000 0.00000 0.00000 


1 0.010 

2 0.026 

3 0.012 
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0.043 
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40 0.042 

41 0.030 
9D 

1 1.2000 1.0000 0.0000 1.0000 - 1.0000 7.5000 

2 1.2720 2.1070 0.0000 1.4500 -1.1500 5.1720 

3 1.4850 3.5250 0.0000 1.8500 -1.8500 3.2760 


7 1.6590 4.1210 0.0000 1.8500 - 1.8500 3.2960 

10 1.5340 3.5040 0.0000 1./300 -1.7500 3.6630 

11 1.8680 5.6450 0.0000 2.1000 -2.1000 2.5510 

SA 

FORHflT=COUPLED 
NHAT= 1 

1 45.5 45.5 45.5 0. 0. 0. ) 

45.5 45.5 45.5 50.0 50.0 50.0 ) 

45.5 45.5 45.5 -50.0 -50.0 -50.0 

0.388E+06 0. 120E+06 0.388E+06 0. 0. 0.134E+06 0. ) 

0. 0. 0.324E+06 0. 0. 0. 0.101E+06 ) 

0.324E+06 0. 0. 0. 0. 0. 0.112E+06 

NMAT= 1 

2 31.3 31.3 31.3 0. 0. 0. > 

31.3 31.3 31.3 29.2 29.2 29.2 > 

31.3 31.3 31.3 -29.2 -29.2 -29.2 

0.565E+06 0. 175E+06 0.565E+06 0. 0. 0.195E+06 0. > 

0. 0. 0.657E+06 0. 0. 0. 0.204E+06 > 

0.657E+06 0. 0. 0. 0. 0. 0.227E+06 

NHAT= 1 

3 50.0 50.0 50.0 0. 0. 0. > 

50.0 50.0 50.0 161. 161. 161. > 

50.0 50.0 50.0 -161. -161. -161. 

0.353E+06 0. 109E+06 0.353E+06 0. 0. 0.122E+06 0. > 

0. 0. 0.352E+05 0. 0. 0. 0.109E+05 > 

0.332E+05 0. 0. 0. 0. 0. 0.122E+05 


i 

NMAT= 2 

154 0. 0. 0. 0. 0. 0. > 

0 . 0 . 0 . 0 . 0 . 0 . > 

0 . 0 . 0 . 0 . 0 . 0 . 

0. 770E+06 0.298E+05 0.192E+06 -10.0 -10.0 0.962E+05 10.0 

10.0 -10.0 0. 157E+06 10.0 10.0 -10.0 0.122E+05 

0.393E+05 -10.0 -10.0 10.0 -10.0 -10.0 0.197E+05 

MMAT= 2 

155 0. 0. 0. 0. 0. 0. ) 

0 . 0 . 0 . 0 . 0 . 0 . > 

0 . 0 . 0 . 0 . 0 . 0 . 

0.385E+06 0. 149E+05 0.962E+05 -10.0 -10.0 0.481E+05 10.0 

10.0 -10.0 0.301E+05 10.0 10.0 -10.0 0.621E+94 

0.200E+O5 -10.0 -10.0 10.0 -10.0 -10.0 0.100E+05 

»1AT= 2 
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136 

0. 

0 . 

0 . 

0 . 

0 . 


0. 0. 

0 . 

0 . 


0. 0. 

0.292E+07 0. 113E+06 

0. 729E+06 -10.0 

to.o 

-10.0 


0.S43E+06 10.0 

0. 136E+06 -10.0 


-10.0 10.0 

COM= 1 




ZERO i: 

189 



ZERO 3: 

189 



ZERO 3: 

237 



SYMMETRY PLANE=2 



RMASS 




ftPATRAN NODAL TEMP SET 

101 

93 

2237.800 



96 

2367.600 



98 

2132.700 



118 

864.600 



121 

964.800 



123 

991.600 



126 

941.800 



4PATRAN 

NODAL TEMP SET 

102 

93 

2223.400 



96 

2427.200 



98 

1310.200 



118 

699.900 



121 

1908.000 



123 

2102.100 



126 

• 

1219.100 



• 

« 

4PATRAN 

NODAL TEMP SET 

240 

146 

585. 25 



160 

385.23 



218 

3279.38 



233 

3279.38 



238 

1824.38 



232 

1824.38 



91 

163.66 



103 

163.66 



117 

806.34 



130 

806.34 



193 

378.29 



214 

878.29 



[XQT ELD 



£23 




Nt1AT= 

l: NSECT= 


i: NNSW= 1 

58 

81 0 

0 *el no 1 

r«AT= 

i: NS£CT= 


l: NNSM= 1 

57 

80 0 

0 #el no 2 

!MAT= 

l: NSECT= 


l: NNSW= 1 

60 

83 0 

0 *ei HQ 3 



0 . 0 . 

) 

0 . 

0 . 

> 

0 . 

0 . 


- 10.0 

0.365E+06 

10.0 > 

10.0 

-10.0 

0.421E+05 ) 

- 10.0 

-10.0 

0.679E+05 
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1 


NMAT= 

l: 

NSECT= 

l: NNSW= 1 


50 

43 

0 

0 tel na 575 


NMAT= 

l: 

NSECT= 

l: NNSM= 1 


48 

42 

0 

0 tel na 576 


:NHAT= 

l: 

NSECT= 

l: NNSW= 1 


47 

E24 

37 

0 

0 tel na 577 


?«1AT= 

l: 

NSECT= 

l: NNSW= 18: 

NREF= 

398 

397 

0 

0 tel na 578 


NMAT= 

i: 

NSECT= 

i: msw= 18 : 

NRtF= 

397 

396 

0 

0 tel no 579 


«1AT= 

i: 

NSECT= 

i: nnsw= is: 

NR£F= 

396 

• 

395 

0 

0 tel na 580 


• 

« 

M1AT= 

l: 

NSECT= 

i: nnsw= 18: 

NREF= 

211 

208 

0 

0 tel na 804 


NMAT= 

i: 

NSECT= 

l: nnsw= is: 

NREF= 

232 

229 

0 

0 tel na 305 


NHAT= 

l: 

NSECT= 

11: NHSW* 28: 

NREF= 

249 

E33 

247 

0 

0 tel na 806 


NMAT= 

l: 

NSECT= 

l: NNSW= 29 


162 

147 

175 

0 tel na 807 


NHAT= 

l: 

NSECT= 

2: NNSW= 30 


161 

145 

163 

0 tel na 808 


NHAT= 

l: 

NSECT= 

3: NHSW= 31 


147 

120 

151 

0 tel na 809 


• 

• 

M1AT= 

2: 

NSECT= 

11: NNSW= 39 


120 

116 

118 

0 tel na 836 


NMAT= 

2: 

NSECT= 

11: NHSM= 39 


95 

92 

93 

0 tel na 837 


NHAT= 

2: 

NS£CT= 

11: NHSW= 39 


73 

E43 

69 

72 

0 tel na 338 


:tWT= 

l: 

NS£CT= 

14: NNSW= 41 


64 

87 

89 

66 tel na 339 


Nf1AT= 

l: 

NSECT= 

8: NNSW= 36 


63 

86 

88 

65 tel na 840 


NMAT= 

l: 

NSECT= 

7: NNSH= 35 


66 

« 

89 

90 

67 tel na 341 


• 

*»• 

NSECT= 

143: NNSW= 39 


189 

210 

207 188 #ei na 1312 


-WAT® 

^ • 
«»• 

NSECT= 

143: NKS«= 29 
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210 231 223 207 *el na 1313 
NMAT= 2: NSECT 3 143: NNSW 3 3? 

231 237 236 22S #el na 1314 
[XQT E 

T= 1.-20. .05. 1.-3. l.-l. 20.. 1.-4, 1.-4, 1.-4 
[XQT EKS 
CXQT TOPO 
[XQT K 
[XQT N 
[XQT AUS 

MAS=*SUM(CEJ1,R«AS) 


[XQT INV 
RESET CON 3 1 
[XQT AUS 

ALPHA: CASE TITL 1 

l’wing tip loading 
SYSVEC: APPL FORC 1 

Js 38; 0. 0. 0.2500E+05 0. 0. 

j- 64: 0. 0. 0.2500E+05 0. 0. 

♦call forces ua (lib 3 <l>: set=( 1>! con 3 U>; jnt : U>> 

[XQT AUS 

ALPHA: CASE TITL 2 

l’aero case 1 
ELDATA: PRES E33 2 1 

S=i: E= 2: 0.495479 0.495479 0.495479 

6=1 : E 3 4: 0.476364 0.476364 0.476364 

6=1; E= 6: 0.459449 0.459449 0.459449 

6=i: E= a: 0.429918 0.429918 0.429918 

Q>i: E= 10: 0.396372 0.396372 0.396372 


ELDATA: PRES E43 2 1 

6=1: E 3 2: 0.372723 0.372723 0.372723 
6=1: E 3 4: 0.367201 0.367201 0.367201 
3=1: E 3 6: 0.460492 0.460492 0.460492 


0.372723 

0.367201 

0.460492 


0 . 

0 . 


5=1: E= 132: 0.448571 0.448571 0.448571 0.448571 

g=l': E= 134: 0.623787 0.623787 0.623787 0.623787 

[XQT EQNF 
RESET SET 3 2 
[XQT DCU 

CHANGE 1 EQNF FORC 2 1 APPL FORC 2 1 

♦call forcesua 1 1 ib = < 1 ) ? set 3 < 2)? con 3 <l>? jnt-il)) 

[XQT AUS 

ALPHA: CASE TITL 3 
i’aero case 2 
3JDATA: PRES E33 3 1 

3=t; E 3 2: 0.664626 0.664626 0.664626 

3=i: E* 4: 0.594172 0.594172 0.594172 

3=15 E= 6 : 0.586386 0.586386 0.586386 

G=i: E= 3: 0.633090 0.633090 0.633090 

-;=i: E= 10: 0.734930 0.754930 0.754930 

ELDATA: PRES E43 3 1 
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■j*i; E= 2: 0.330738E-01 0.530738E-01 0.5507381- -01 0.550738E-01 
j»l: E= 4: -0.919826E-01 -0.9198266-01 -0. 919826E-01 -0. 919826E-01 


3=1: E= 6: 

• 

0. 129170 

0.129170 

0.129170 

0. 129170 

• 

3=1: E= 130: 

0.869722 

0.869722 

0.369722 

0.869722 

3=1: E= 132: 

0.983188 

0.983188 

0.985188 

0.983188 

3*1: E= 134: 

0.789944 

0.789944 

0.789944 

0.789944 

CXQT EQNF 





RESET SET= 

3 




CXQT XU 





CHANGE 1 EQNF 

FORC 3 1 APPL FORC 

3 1 



♦call forces ue (lib*(l)i set=< 31! con— < 1 ) » j nt - < 1 > ) 


CXQT AUS 

ALPHA: CASE TITL 4 

l'aero case 3 
eldata: PRES E33 4 1 
3=1: E= 21 0.336062 

0.336062 

0.336062 


8-1 : E= 4: 0.276674 

0.276674 

0.276674 


3=1 : E= 6: 0.249004 

0.249004 

0.249004 


3=1 : E= 8: 0.233692 

0.233692 

0.233692 


8=1 : E= 10: 0.245430 

0.245430 

0.243430 


ELDATA: PRES E43 4 1 
8=1: E= 2: 0.139290 

0.139290 

0.139290 

0.139290 

8=1: E= 4: 0.969320E-01 0.969320E-01 

0.969320L-01 0.969320E-01 

3=1: E= 6: 0.960526E-01 0.960526E-01 0.960526F-01 0.960526E-01 

• 

t 

8=1 : E= 130: 0.307940 

0.307940 

0.307940 

0.307940 

G=i: E= 132: 0.312033 

0.312035 

0.312035 

0.312033 

3=1: E= 134: 0.412321 

0.412321 

0.412321 

0.412321 

CXQT EQNF 
RESET SET= 4 
CXQT DCU 

CHANGE 1 EQNF FORC 4 

1 APPL FORC 

4 1 


♦call forcesu* <lib=<l) ! 

set=< 4); con 

=<l>; jnt* 

<1>> 

CXQT AUS 

ALPM: CASE TITL 5 

l’aero case 4 
ELDATA: PRES E33 5 1 
3=1: E= 2: 1.51809 

1.51809 

1.51309 


3=1 : E= 4: 1.44913 

1.44913 

1.44913 


3=1: E= 6: 1.40116 

1.40116 

1.40116 


3=1: E= 3: 1.29837 

1.29837 

1.29337 


3=1: E= 10: 1.13715 

1. 13715 

1.13715 


ELDATA: PRES E43 5 1 
3=1: E= 2: 0.983338 

0.983358 

0.983358 

0.983358 

3=1: E= 4: 0.971026 

0.971026 

0.971026 

0.971026 

1=1: E= 6: 1.37643 

1.37643 

1.37643 

1.37643 
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3=i : e= 130: 

1.62700 

1.62700 

1.62700 

1.62700 

3=t: E= 132: 

1.38076 

1.38076 

1.38076 

1.38076 

3=1 : E= 134: 
[XQT EQNF 

2.01220 

2.01220 

2.01220 

2.01220 

RESET SET= 
[XQT DCU 

5 




CHANGE 1 EGWF FORC 5 1 APPL FORC 

5 1 


♦call forcesua (Iib=<l)i 

set=( 5); 

can=(l>; j nt- <1>) 


[XQT AUS 





ALPHA: CASE TITL 6 

l’aero case 5 




eldata: PRES 

E33 6 1 




o=i : E= 2: 

1.63482 

1.63482 

1.63482 


G=i: E= 4: 

1.36996 

1.36996 

1.36996 


s=i: E= 6: 

1.21555 

1.21555 

1.21555 


G=i: E= 8: 

1.09577 

1.09577 

1.09577 


G»i: E= 10: 

1.07144 

1.07144 

1.07144 


eldata: PRES 

E43 6 1 




g=i: e= 2: 

0.609049 

0.609049 

0.609049 

0.609049 

G=i: E= 4: 

0.405767 

0.405767 

0.405767 

0.405767 

3=1 : E= 6: 
% 

0.469569 

0.469569 

0.469569 

0.469569 

m 

• 

G=i: E= 130: 

1.47321 

1.47321 

1.47321 

1.47321 

G=i: E= 132: 

1.37584 

1.37584 

1.37584 

1.37584 

G*i: E= 134: 

2.15222 

2.15222 

2.15222 

2.15222 

[XQT EQNF 
RESET SEJ- 
[XQT DCU 

6 




CHANGE 1 EQNF FORC 6 

1 APPL FORC 

6 1 



•call forces ui < 1 ib— < 1 > » set=< 6); con=<l)i jnt r <l)) 
[XQT AUS 

ALPHA: CASE TITL 50 
l’l.Og inertial load 
[XQT AUS 
R=RIGID( 225) 

DEFINE F=1 R AUS 1 1 3 3 
NX2 a PRQDUCT (MAS.F) 

APPL FORC 50 1 =UNION( -1.00 MX2) 

•call forcesua <lib=U); set=< 50); con=(l>; jnt= <225) ) 
[XQT AUS 

ALPHA*. CASE TITL 60 

l’2.5g inertial load 
CXQT AUS 
R=RIGID< 225) 

DEFINE F=1 R AUS 1 1 3 3 
1X2=PRQDUCT (MASiF) 

VPL FORC 60 1 =UNION( -2.50 MX2) 

•call forces ua ( 1 ib = * 1 ) > set=( 30); ~cn=M ) '■ 1 ) 

.XQT AUS 

iLPHA: CASE TITL 51 
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1 1 load case 1 
CXQT aus 

DEFINE AA=1 APPL F0RC 2 1 

DEFINE BB=1 APPL F0RC 50 1 

APPL FORC 51 1 =SUf1( 1.0000 AA, 0.7311 BO) 

♦call forcesua (lib=(l>! set=( 51); con=(l); jnt-(l)) 
CXQT AUS 

ALPHA: CASE TITL 52 
IMoad case 3 
CXQT aus 

DEFINE AA=1 APPL FORC 31 

DEFINE BB=1 APPL FORC 50 1 

APPL FORC 52 1 =SUM( 1.0000 AA, 0.9708 B6) 

♦call forces us (lib=<l)i set=< 52); con=<l>; j nt ~ < 1 > ) 
CXQT AUS 

alpha: CASE TITL 53 

l’load case 3 
CXQT aus 

DEFINE AA=1 APPL FORC 4 1 

DEFIfE BB=1 APPL FORC 50 1 

APPL FORC 53 1 =SUM< 1.0000 AA, 0.4928 SB) 

♦call forcesua (lib=(l); set=< 53); con=(l>: jnt-(l)) 
CXQT AUS 

alpha: CASE TITL 61 

IMoad case 4 
CXQT aus 

DEFINE AA=1 APPL FORC 51 

DEFINE 8B=1 APPL FORC 60 1 

APPL FORC 61 1 =SUN( 1.0000 AA, 0.9220 BB) 

♦call forcesua (lib=<l>! set=< 61); con=(l); j nt = < 1 > ) 
CXQT AUS 

ALPHA: CASE TITL 62 
IMoad case 5 
CXQT aus 

DEFINE AA=1 APPL FORC 61 

DEFIIC BB=1 APPL FORC 60 1 

APPL FORC 62 1 =SUM( 1.0000 AA, 0.9827 BB) 

♦call forcesua (lib=(l); set=< 62); con=<l); j nt = < 1 > ) 


CXQT SSOL 



RESET CON= 
CXQT SSOL 

1 SET= 

51 

RESET CON= 
CXQT SSOL 

1 SET= 

52 

RESET CON= 
CXQT SSOL 

1 SET= 

53 

RESET CON= 
CXQT SSOL 

1 SET= 

61 

-ESET CON= 
XQT EXIT 

1 SET= 

62 
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Appendix E 


Iterative Procedure for Balancing Airplane 
Load Case No. 1 

The process used to balance the loads on the airplane for 
trimmed flight is described in this Appendix. First, an initial 
structural resizing analysis was conducted to establish the total 
weight (structural and non-structural) of the airplane and the lift 
component of the aerodynamic pressures. Then, three or more 
iterations were performed to adjust the aerodynamic lift and 
balance load for a trimmed airplane. 

Initial Execution 

This initial COMET execution is for determining the lift 
component of the aerodynamic pressures and the magnitude of the 
total weight so that factors can be determined for matching these 
values to the design loads. The lift component and total weight are 
determined during the execution by a procedure called "forcesum". 
This procedure determines the summation of forces and moments 
about a specified reference point for a specified set of loads (the 
reference point is node 225 for the present study). These 
resultants, for the symmetrical half of the airplane, are written to 
an output (or log) file that is created during program execution and 
are as follows: 

Aerodynamic Lift (load case 2) 202,500 lbs. 

Inertial Force (load case 50) 381,400 lbs. 

The weight fraction used for the initial resizing exercise is 
determined as follows using the flight design gross weight and 
maximum gross weight for a total airplane given in Tables A-1 and 
A-2 of Appendix A, respectively: 

Design GW/Maximum GW = 494,000/675,000 = 0.7311. 

Multiplying the inertial force by the weight factor and adding the 
aerodynamic lift results in an imbalance of 152,640 lbs. (load set 
51) in the downward direction; which requires an increase in the 
aerodynamic lift to correct. However, since the non-structural 
weights of Appendix B were used in formulating the finite element 
model, it was decided that the corresponding maximum gross weight 
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of 713,696 lbs. would be the baseline weight for the present study. 
Therefore, a revised design gross weight for succeeding iterations 
i s 

(.731 1)(713,696) = 521,783 lbs. 

First Iteration 

This COMET execution is performed to obtain a balance of the 
inertial force and the aerodynamic lift component. Using the 
inertial force (which is based on a 1.0 g maneuver) of the initial 
execution, the calculated maximum gross weight is (2)(381,400) = 
762,800 lbs. The new weight fraction for this iteration is 

Design GW/Maximum GW = 521,783/762,800 = 0.6840. 

The aerodynamic lift, as given by the initial execution, is 
(2)(202,500) = 405,000 lbs. Therefore, the aerodynamic fraction 
required to balance the forces is 

Design GW/ Aerodynamic Lift = 521,783/405,000 = 1.2884. 

Second Iteration 

The unbalanced moment about the Y-axis at node 225 due to the 
offset of the center of gravity and the center of pressure is given by 
the first iteration (load set 51) as 21,100,000 in.-lbs. A balance 
load is applied at node 391 (an arbitrary selection) which is 1531.46 
inches aft of node 225. The balance load is 

21,100,000/1531.46 « 13,778 lbs. 

To account for this additional vertical force, the aerodynamic 
fraction is modified as follows: 

[521,783 - (2)(1 3,778)]/405,000 = 1.2203. 

Third Iteration 

The reduction of the aerodynamic lift by the balance load 
causes a secondary change in the moment imbalance between the lift 
and the inertial loads. The unbalanced moment from the second 
iteration is -278,500 in.-lbs which requires a reduction in the 
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balance load by 182 lbs (278,500/1531.46) to 13,596 lbs (13,778 - 
182) and another modification of the aerodynamic fraction as 
follows: 

[521,783 - (2)(1 3,596)] / 405,000 = 1.2212. 

No further iterations were performed. 
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